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ABSTRACT 


In  this  paper  we  consider  the  laminar  boundary  layer  which  forms  on  a  sharp  cone  in  a 
supersonic  freestream,  where  lateral  curvature  plays  a  key  role  in  the  physics  of  the  problem. 

This  flow  is  then  analysed  from  the  point  of  view  of  linear,  temporal,  inviscid  stability. 
Indeed,  the  basic,  non-axisymmetric  disturbance  equations  are  derived  for  general  flows  of 
This  class,  and  a  so  called  “triply  generalised”  inflexion  condition  is  found  for  the  existence  of 
i  'J,subsomc’r  neutral  modes  of  instability.  This  condition  is  analogous  to  the  well-known  gen¬ 
eralised  inflexion  condition  found  in  planar  flows,  although  in  the  present  case  the  condition 
depends  on  both  axial  and  aximuthal  wavenumbers. 

Extensive  numerical  results  are  presented  for  the  stability  problem  at  a  freestream  Mach 
number  of  3.8,  for  a  range  of  streamwise  locations.  These  results  reveal  that  a  new  mode  of 
instability  may  occur,  peculiar  to  flows  of  this  type  involving  lateral  curvature. 

Additionally,  asymptotic  analyses  valid  close  to  the  tip  of  the  cone  /  far  downstream  of 
the  cone  are  presented,  and  these  give  a  partial  (asymptotic)  description  of  t-Hp  additional 
mode  of  instability.  >  /  V 
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It  is  well  known  that  unlike  many  of  their  incompressible  counterparts, 
supersonic  boundary  layers  arc  generally  susceptible  to  inviscid  forms  of 
instability;  this  comment,  in  particular,  applies  to  Blasius-type  boundary 
layers.  The  paper  which  first  tackled  the  problem  of  stability  of 
supersonic  boundary  layers,  from  any  kind  of  rigorous  mathematical 
standpoint  was  that  of  Lees  and  Lin  (1946),  who  presented  some  extremely 
important  results  on  the  problem,  including  the  importance  of  the  quantity 

— *  |p*  — *1  (where  u*  denotes  velocity  tangential  to  the  surface,  v* 
dy  1  dy  > 

the  normal  to  the  surface,  and  p*  the  fluid  density),  which  plays  the 
same  role  as  that  of  an  inflexion  point  in  incompressible  flows.  In 
particular  at  the  point  where  the  above  expression  is  zero  (y*  =  y0*.  say), 
then  a  neutral  mode  exists  with  wavespeed  u*(y0*)  (although  see  the 
restriction  (4.21)  later  in  this  paper).  Lees  and  Lin  also  showed  how 
viscosity  played  a  similar  role  at  a  critical  layer  as  in  analogous 
incompressible  flow  situations.  This  two-dimensional  work  was  later 
extended  to  three -dimens ions  by  Reshotko  (1962). 

The  development  of  digital  computers  enabled  accurate  numerical 
computations  of  the  problem  to  be  made.  One  of  the  earliest  was  that  of 
Brown  (1962),  and  this  was  followed  by  a  number  of  studies  by  Mack 
(1963,  1964,  1965a, b,  1969,  1984,  1987a),  both  of  an  inviscid  (large 
Reynolds  number)  and  viscous  (finite  Reynolds  number)  nature.  These 
computations  revealed  a  further  distinction  between  supersonic  and 
incompressible  flows,  namely  that  in  the  latter  case  an  infinite  sequence  of 
modes  is  possible.  These  are  important,  because  generally  it  is  one  of 
these  "higher"  modes  that  exhibits  the  largest  growth  rate  according  to 
inviscid,  linear  stability  theory. 

In  spite  of  their  great  practical  importance,  flows  of  this  type,  but 
involving  lateral  curvature  have  received  scanty  attention  over  the  years. 
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Duck  and  Hall  (1989a)  showed  how  in  supersonic  flows,  curvature  interacting 
with  viscosity  could  provoke  additional  instabilities  (ax i symme t r i c  in 
form),  provided  the  body  radius  was  below  some  critical  value.  Duck  and 
Hall  (1989b)  then  went  on  to  show  how  a  similar  effect  occured  with 
non-axi symme t ri c  modes  (which,  in  fact,  turn  out  to  be  generally  more 
unstable  than  corresponding  axisymmetric  modes). 

Recently,  Duck  (1989)  (hereafter  referred  to  as  I)  has  shown  the 
influence  of  curvature  on  the  axisymmetric  inviscid  stability  of  supersonic 
boundary  layers,  in  particular  those  that  form  on  thin  straight  circular 
cylinders.  It  was  found  that  curvature  has  a  profoundly  stabilising  effect 
on  these  modes.  Interestingly,  it  was  shown  how  curvature  altered  the 
generalised  inflexion  condition  described  above,  and  a  modified  (or  "doubly 
generalised")  inflexion  condition,  involving  the  radius  of  curvature  was 
derived. 

In  the  present  paper  we  extend  the  work  of  I  to  non-axi symmet ric- 
disturbances  which,  indeed,  turn  out  to  be  more  important  than  axisymmetric 
disturbances  considered  previously.  Further,  rather  than  studying/applying 
our  techniques  to  the  thin  straight  circular  cylinder,  we  consider  a 
somewhat  more  practical  configuration,  namely  that  of  a  sharp  cope. 

To  (considerably)  facilitate  already  lengthy  and  complex  computations, 
following  I  we  assume  that  the  lip  of  the  cone  has  a  finite  radius,  and 
that  the  associated  boundary  layer  on  the  cone  surface  is  planar  at  the  cone 
tip.  (We  expect  that  far  downstream  the  flow  will  become  planar  once  again, 
having  utilised  the  Mangier,  1948,  transformation.)  Following  previous 
studies  in  this  area,  we  also  ignore  the  presence  of  any  shocks. 

Mack  (1987b)  has  performed  some  computations  for  the  stability  of  the 
flow  over  a  cone  in  supersonic  flow,  at  finite  Reynolds  numbers,  but  found 
little  difference  with  corresponding  planar  results.  Here,  we  deliberately 
allow  curvature  to  occur  throughout  the  study,  both  in  the  equations 
governing  the  basic  flow,  and  in  the  disturbance  equations. 
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A  few  preliminary  numerical  results  of  this  study  have  been  presented  by 
Duck  (1990).  Here,  we  present  additional  numerical  results,  and  also  a 
number  of  important  asymptotic  results  related  to  these  calculations. 


2.  Equations  of  motion  and  state 

The  general  layout  of  the  problem  is  shown  in  Fig.  1.  The  z*  axis 
lies  along  the  cone  axis,  r*  denotes  the  radial  coordinate,  and  0  the 
azimuthal  coordinate.  The  velocity  vector  v*  has  components  vj*.  V2* 
and  V3*  in  the  r*.  9  and  z*  directions  respectively.  Although  we 

shall  be  concerned  with  a  basic  flow  which  is  independent  of  0  (and  with 
V2*  =  0),  when  we  go  on  to  consider  the  stability  of  the  flow,  we  shall  be 
concerned  with  non- ax i symmetric  disturbances. 

In  the  cylindrical  polar  coordinate  system  as  defined  above,  the  full 
equations  of  continuity,  momentum  and  energy  take  on  the  following  forms 


(Thompson  1972) 


dp*  +  a 


.  -  ,  .  (P*  vl*)  +  K  (P*  v2*> 

dt  dr  r  30 

♦  * 

P  vj  _ 


/  *  r  1 1  _  n 

+  T%(P  v3  )  +  v~  -  Ul 
dz  r 


(2.1) 


*  f  Dvl*  (v2*) 


Dt 


21  =  + 

3r*  3r 

+  I* 

r  30  x  dz 

+  h*7.-.  M  , 


(2.2) 


*  r  Dv^  V1*V2*  1  =  .  dp*+ 

H  L  DT  r  J  7  30  3r 


+  I_  9  -99  .  9  gfl L 
7  30  3z 

+  2  lr*0 


(2.3) 
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=  .  V  +  Lhlil  +  '  3 

r\  .  *  *  i  *  *  r\ 


di  dr 


r  90 


+  ;>  Y  *  *  Y  *  * 

+  ?  kl^L  +  i^_£_ 
9z  r 


(2.4) 


*  D  - 

p  ^(cp 


t*)  -  =  r*  +  ip  fL^  (  kV  ^  ] 

+  —>  L  (  K*  | 

(TV  90  l  90  J 
+  a  (  k*9t*  ] 

+  ^  l  al*  J' 


(2.5) 


Here  p*  is  the  density  of  the  fluid,  p*  the  pressure,  Cp  the 
specific  heat  (at  constant  pressure),  and  K*  the  coefficient  of  heat 


conduction.  The  Eulerian  operator  is  defined  as 


D  9  *  9  vi*  9  •  9 

Dt*  9t*  9r  r  90  9z 


(2.5) 


and  the  viscous  stress  components,  assuming  Newtonian  flow  are  defined 


to  be 


IrV  -  2n*  X*  V .  v*  , 

dr 


*••**[;»  H  ♦  4‘ *-*• 


XZV  =  2p*  ^  +  X*  V.v*  . 

9z 


^  v  *  *  T  1  9vi*  *  9  [v?* 

Ir-e  =  l0r*  =  V  3?V 

v  *  v  *  r  av  i  9vi*  ] 

I92-  =  Iz*p  =  p  [  *-7t  Te"  J  ' 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


v**_v**  *T  ay  avj*  ■ 

iz  ,  -  2-r  z  -  p  "aT*  -  “5^  J 


(2.11) 


The  dispersion  function  f*  in  (2,5)  is  given  by 
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r*  =  2n 


D?V  +  +  D  2/z* 


+  2  D?*q  +  2  D +  2  D2*r* 


+  (X*  -  2  ji*)  (V.v*)2. 

I 


(2.12) 


Here  the  'D'  terms  are  the  components  of  the  rate-of-deformat ion  tensor, 
for  example 


(2.13) 


3v  *  i  3v_* 

- +  — gp  - - 

dz  T7  90 


(2.14) 


The  coefficients  p*  and  X*  above  denote  the  first  coefficient  of 
viscosity  and  bulk  viscosity  respectively  (which  are  expected  to  be 
functions  of  temperature). 

The  equation  of  state  is  taken  to  be  that  pertaining  to  a  perfect  gas. 


i .  e . 


p*  =  p*  R*  T*.  (2.15) 

where  R*  is  the  gas  constant. 

With  reference  to  Fig.  1,  the  surface  of  the  cone  is  taken  to  lie 
along  r*  =  a*  +  X j  z*.  z*  >  0,  (later,  important  assumptions  regarding 
the  size  of  the  slope  parameter  will  be  made),  and  so  on  this  surface  we 
requi re 

vj*  =  V2*  =  V3*  =  0.  (2.16) 

If  the  surface  of  the  cone  is  insulated  then  the  following  additional 
boundary  condition  must  be  imposed 

=  0  ,  (2.17) 

3n 

(where  n*  denotes  an  outwards  normal  to  the  wall).  In  the  case  of  heated/ 
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cooled  walls,  then  t he  condi lion 

T*  =  T*  (2.18) 

w 

must  be  imposed  at  ihc  surface. 

Conditions  remain  to  be  specified  at  z*  =  0;  for  this  we  follow  I 
(precisely)  by  assuming  that  the  boundary  layer  at  this  location  has  zero 
thickness,  enabling  planar  conditions  to  be  imposed  at  this  position. 
Assuming  the  cone  to  be  slender,  then  the  far-field  conditions  are  taken  to 
be  uni  form,  with 


< 

* 

II 

< 

* 

II 

0 

(2.19) 

*  If  ♦ 

V3  =  LL 

(2.20) 

T*  =  T** . 

(2.21) 

(Indeed,  in  the  case  of  non-slender  cones,  if  the  shock  wave  is  attached, 
downstream  of  the  shock  the  external  flow  is  i rrotat ional ,  and  there 
is  also  a  constant  slip  velocity  at  the  surface  of  the  cone.)  We  next  go  on 
to  derive  the  basic  (boundary-layer)  flow  on  the  surface  of  the  cone, 
assuming  curvature  plays  a  key  role  in  the  physics  of  the  problem. 


We  define  our  Reynolds  number  on  the  tip  radius  of  the  cone, 


fol lows 


(3.1) 


and  this  will  be  taken  to  be  large  throughout  this  paper.  As  noted 
previously,  the  basic  flow  is  taken  to  be  independent  of  0,  and  has  no 
azimuthal  velocity  component  (i.e.  V2*  =  0). 

It  is  now  convenient  to  introduce  non-dimensional  parameters  as  follows 
( v j ,  V3,  r,  z,  T,  p,  p)  = 


r»  ♦  *  *  *  T**  *  * 

Re  V]  V3  r  z  T  p  |i 

U„*  ’  LL*  '  a7  ’  Rea*  '  T„*  ’  pM*  ’ 
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A  key  element  of  this  paper  (as  in  1)  is  the  inclusion  of  curvature 


terms  in  the  governing  equations  to  leading  order.  If  this  is  to  be  the 
case,  then  we  must  have  that 


X]  =  Re' *  X 

where  X  =  0(1),  (3.3) 

implying  a  slender  cone. 

The  leading  order  governing  equations  may  then  be  written  (assuming 
Re  — >  ~) 


If  ft)* 

v  1  9  1 

7T  +  37 

II 

O 

(3.4) 

5f-» 

(3.5) 

T  3 
"  7  37 

(3.6) 

*1  37  *  *3  H  =  kT<y-D  M.2  [|fl]2 


where  the  result 


has  been  used,  o  is  the  Prandtl  number,  namely 

n  —  _ *_ 


(3.7) 


(3.8) 


(3.9) 


(which  is  taken  to  be  a  constant  in  this  paper),  y  denoting  the  ratio  of 
specific  heats,  and  is  the  f reest ream  Mach  number,  namely 

M»  =  U„*  /  (y  R*  TM*)i  .  (3.10) 


The  boundary  conditions  to  be  applied  to  this  system  are 
vj  =  V3  =  0  on  r  =  1  +  Xz  , 

V3  — *  1 ,  T  — ♦  1  as  r  — *  »  .  (3.11) 

together  with  a  wall  temperature  condition;  in  the  case  of  insulated  wails, 
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(lo  leading  order) 

on  r  =  1  +  Xz  .  (3.12) 

whilst  in  the  ease  of  healcd/coolcd  surfaces 

T=TW  on  r  =  1  +  Xz  .  (3.13) 

To  close  the  problem,  a  viscosity/temperature  law  is  required.  For  the 
purposes  of  this  paper  we  take  the  linear  Chapman  law  (Stewartson  1964), 
namely 

p  =  CT  (3.14) 

where  C  is  taken  to  be  constant  (although  as  noted  in  I,  there  would  be 
no  conceptual  difficulties  in  taking  more  complex  variations  of  viscosity 
with  temperature). 

If  we  write 


vi  =  C  v,  § 

-1 

z  =  C  z 

X  =  CX  ,  ’  (3.15) 


whilst  retaining  other  terms,  then  the  system  (3.4)  -  (3.7)  becomes 


a 

“Jr 


^  +  rf +  r? 


=  o  , 


Bp  _ 


J7 

9v3 

JT  + 


=  0  , 


3v3  T  a  _  dv3' 

v3  ar  =  7  37  [rT  arj  • 


ay3i 


_  3T  3T  t2  n  u  2  P*3 

V1  JT  +  V3  31  =  T  (Y'U  37- 


T  a  ("rT  3T  " 

+  7  37  a  37 

. 


(3.16) 

(3.17) 

(3.18) 


(3.19) 


whilst  the  wall  boundary  conditions  are  to  be  applied  on 


r  =  1  +  1  2. 
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As  described  previously  we  assume  planar  conditions  prevail  at  z  =  0, 
where  the  boundary  layer  is  taken  to  have  zero  thickness.  The  problem  is 
thus  singular  at  z  =  0,  and  consequently  sealed  variables  must  be 
introduced  in  order  to  (numerically)  solve  the  system  (3.16)  -  (3.19) 
accurately.  Specifically  we  write 


vi  =  C*1  vj  (n,o. 
v3  ='v3  (H.O. 

T  =  T  Ol.O,  (3.20) 

where 

C  =  2*  .  (3.21) 

and  1)  =  (r  -  1  -  XC2)/C-  (3.22) 

The  quantities  with  a  hat  are  then  expected  to  behave  regularly  as 
C  — *  0.  and  indeed  approach  the  planar  solution.  Equations  (3.16)  - 
(3.19)  now  take  the  following  form 


C  3 

IZC 


(3.23) 


(3.24) 


(3.25) 


*  9T  t  Cv 3  3T  “hr  Til  9T 
V1  +  ^  3C  *  v3  [K  *i\  ^ 


=  t2(y-i)m«2 


2 

Ov-3  T9 

F  rT  9T] 

+  7 

[o' 

(3.26) 
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Here  of  cours< 

r  =  I  +  XC2  +  Cn  .  (3.27) 

The  boundary  conditions  in  terms  of  these  variables  arc 
v  1  =  v3  =  0  on  T]  =  0 

V3  — »  1,  T  — »  1  as  t)  — *  «  .  (3.28) 

In  the  case  of  insulated  walls,  the  additional  surface  condition  is 

gl  =  0  on  t)  =  0  .  (3.29) 

whilst  for  heatcd/cooled  walls, 

T  =  Tw  on  T)  =  0  .  (3.30) 

Setting  £  =  0  reduces  the  system  to  an  ordinary  differential  system 
(corresponding  to  the  planar  case)  as  in  I  in  the  same  limit  (indeed, 
setting  X  =  0  reduces  (3.23)  -  (3.26)  to  the  corresponding  system 
considered  in  1).  The  solution  of  the  ordinary  differential  system  at 
£  =  0  then  provides  the  initial  conditions  for  a  (straightforward) 
Crank-Nicolson  scheme  in  £,  identical  to  that  used  in  I.  For  the 
purposes  of  this  paper,  we  shall  focus  our  attention  on  the  insulated  wall 
case  (although  the  heated/cooled  case  may  be  treated  in  exactly  the  same 
manner).  Distributions  of  wall  temperature  with  axial  coordinate  £(=  ii) 
are  shown  in  Fig.  2a  (M*,  =  2.8)  and  Fig.  2b  (M„  =  3.8),  and  the 
corresponding  distributions  of  wall  shear  V3fj i ti  =  0  are  shown  in  Fig.  3a 
(Mo.  =  2.8)  and  Fig.  3b  (M^  =  3.8). 

In  all  cases,  these  distributions  are  quite  different  to  the 
corresponding  X  =  0  distributions  (i.e.  the  distributions  on  a  straight 
circular  cylinder)  as  found  in  I,  where  at  both  =  2.8  and  3.8 

a  monotonic  decrease  in  values  was  found.  It  is  also  quite  clear  that  the 
results  evolve  from  the  planar  case,  to  the  far  downstream  limit,  as 
predicted  by  the  Mangier  Transformation  (Mangier  1946,  Stewartson  1964), 
name  1  y 
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V^r  v3ti  n  =  o.  C  =  0 


T 


C  — *  " 


C  =  0  . 


(3.31) 


In  the  following  section  we  go  on  to  investigate  the  stability  of 
flows  of  this  class,  subject  to  small  amplitude  inviscid  disturbances. 


_ Inviscid  s lability  of  the  flow 

U _ Disturbance  equations 

In  this  section  we  derive  the  disturbance  equations  relevant  to  small 
amplitude  disturbances  in  any  supersonic  axisymmetric  boundary  layer  type 
flow.  In  1,  just  axisymmetric  disturbances  were  considered;  here  we 
consider  the  more  general  case  of  non-axi symmet ric  perturbations  of  the 
f  1  ow . 

We  consider  disturbances  whose  wavelength  in  the  axial  direction 
(1/a)  is  comparable  to  the  (tip)  radius  of  the  cone.  Specifically,  at  a 
fixed  z  stat ion  we  wri te 

vj*  =  8  ft  IL*  V]  (r)  E  +  0(52)  , 

v2*  =  5  U„*  v2  (r)  E  +  0(5-)  , 

V3*  =  IL*  [W0  (r)  -t  ?  v3  (r)  E]  +  0(62)  , 

T*  =  T^*  [T0  (r)  +  5  T  (r)  E]  +  0(52)  , 

P*  =  P»*  ^  (r)  +  5p  (r)  E]  +  0(52)  . 

P  =  P-*  R*  T„*  [1  +  6p  (r)  E]  +  0(52)  .  (4.1) 

where  E  =  exp  [i  a  (z  -  ct)  +  in  0}  ,  (4.2) 


and  5  is  the  scale  of  the  disturbance  (taken  to  be  dimini  shingly  small), 
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whilst 


t  =  (U«*/a*)  t*  . 
z  =  z*/a* 


wo(0  =  v3  (r,z)  , 
T0(r)  =  T  ( r , z)  , 


(4.3) 


where  V3  and  T  are  found  from  the  computations  in  the  previous  section. 

Substituting  (4.1)  into  (2.1)  -  (2.4)  and  taking  the  0(5)  terms  of 
leading  order  in  Re  yields  the  following  linear  system 


.  ~  i V3  1  Vt ,  v? 

lcp  +  “To  +  rTo  Vl  +  To  +  tn  r~To 


+  i  W0p  -  v  1  =  0 


(4.4) 


ic  -  i  Wn  vi  ~  Wnr  -  ip 
—  V3  +  -2 — a.  +  vi  -Hi.  =  — £.  , 


—  V3  + 

To 


_u — 1  +  v.  *  — * 

To  1  To  yMi 


(4.5) 


i£  vi  +  1  ^2  -  '  »"P 

To  1  To  ' 


(4.6) 


i  a2  c  vi  i  a2  w0 
To  +  To 


=  -  £v 


(4.7) 


To  ['  ^  T  +  iW0  T  +  vj  Tor]  +  Ill  (i 


icp  -  iW0p)  =  0  . 

(4.8) 


The  perturbation  equation  of  state  is 


p  =  To  p  + 


(4.9) 


After  some  algebra,  this  system  may  be  reduced  to  the  following  two 


equa  t i ons 


1  -  Wor  - 

vlr  ♦  7  V1  '  TVoVl 


-  M? 

”  yM2«(W 


Tof  1  +  -^L  ■  M20,(W0-c)2  .  (4. 

o-c)  (  L  a2r2 


■—  (wo-c)  M  =  -  775-  Pr 
To  yM*-« 


(4.11) 
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a  result  also  to  be  found  in  t he  work  on  the  stability  of  jets  by  Michalkc 
(1971).  Writing 


V!  =  C«P.  4  =  a/C  . 


(4.12) 


gives 


C  wnt1  <j> 

-  <P  -  -QD-2 

i+xC2+Cn  w0-c 


- (  To  [l  +  -» — ■  M-2(Vc)2  • 

M2„(W0-C)  [  L  a2  (l+xc2  +  Cn)2  J 


together  wi th 


ia2  (W0-c)  2_  =  '  £lL_  . 

To  yM2* 

These  equations  may  be  combined  to  eliminate  p  to  give 


(4.13) 


(4.14) 


d  (w°‘c)  C  9n  +  'f7Xc2'+'Cn  ^  '  WoT1  9 

dT)  To  1  +  — n-2-^L - ,  ]  -  M2  (W0-, 

a2  (i+x  C2  +  Cn)2  J  °°  0 


-  *1  (V^)2 


-  73  (W0-c)  <P  ■ 

Alternatively  cp  may  be  eliminated  to  give 

~p{(Vc)2  -  TO  [  I  »  a2  (,  ^  ]  } 


(4.15) 


W0T1  .  -  Eft— la — 

0T1  1  +  AC2  +  Cn  a2  (W0-c) 


-  (w0-c)  —  Eu — 

dq  (a2  (W0-c) 


(4.16) 


These  equations  reduce  naturally  to  the  system  considered  in  I  by 
setting  n  =  0  (i.e.  axisymmetric  disturbances)  and  X  =  0  (zero  cone 
angle) . 


The  boundary  conditions  to  be  applied  to  the  above  systems  are  that 
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ip  =  Pr|  =  0  on  T|  =  0  .  (4.17) 

togc  the  r  with 

» -  •  *- {*..♦1  <">  *  Vn  <">}  •  (4I8) 

p  -  T  8^1  1  *n  (I»  .  (4.19) 

[1  -  M„2  ( 1  - c ) 2 ] i 

where  r|  =  ±  a[l-M„,-  ( 1  - c ) 2] ^  +  X(,  +  Tij  ;  (4.20) 

here  the  appropriate  sign  is  chosen  to  ensure  that  the  real  part  of  q 
is  positive  as  q  — »  »>,  to  ensure  boundedness.  Equations  (4.15)  together 
with  (4.17)  and  (4.18)  (or  (4.16)  together  with  (4.17)  and  (4.19),  or 
(4.13),  (4.14)  together  with  (4.17),  (4.18),  (4.19))  constitute  an 
eigenvalue  problem;  here,  we  consider  just  the  temporal  problem,  i.e.  given 
n  and  a  (real),  to  find  c  (generally  complex).  The  problem  was  treated 
numerically  using  a  Runge-Kuita  method  based  on  (4.13)  and  (4.14),  shooting 
inwards  from  some  suitably  large  value  of  tj .  with  contour  indentation  in 
the  manner  described  in  I. 

4.2  The.  :triply  generalised  inflexion  CQndiim". 

In  I,  a  condition  for  the  existence  of  so  called  inviscid,  neutral, 
axisymmetric  mode,  i.e.  those  for  which 

1-1/M*  <  c  <  1  +  1/M„,  (4.21) 

implying  disturbances  that  decay  as  q  — *  «>,  was  derived.  Here,  we  go  on 
to  derive  the  corresponding  condition  for  non-axi symme t r i c  modes. 

If  we  multiply  (4.15)  by  q>*  (where  an  asterisk  denotes  a  complex 
conjugate),  and  subtract  from  the  resulting  equation  its  complex  conjugate, 
then  the  following  equation  is  obtained 
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<p*  d 
W0-c  Hq 


(VO  cpn  +  7  <p>  -  V  v 
x 


<P  d  (W0-c*)  (tpr,*  +  -p  <p*)  -  W011  ip* 
w0-c*  dn  Xs 

where  we  have  written 


X  =  To 


+  ^71  }  '  m2“  (wo-02  . 


and  r  =  ( 1  +  +  ^ri)  /C  • 

After  some  algebra,  (4.22)  can  be  written  as 


* 


<P 


d 

dri 


*  1  * 

<P  T1  +  7  <P 

- f - 

X 


(4.22) 


(4.23) 

(4.24) 


=  r  <p<p* 


1  d 
Vc 


Writing  c  =  cr  +  i  c j  , 

then  the  neutral  state  corresponds  to  Cj 

In  this  limit,  (4.25)  may  be  written 


0. 


(4.25) 

(4.26) 


=  2i  r2  1  <p  1 2c j  d 
lW0-c  12  dq 

Following  the  arguments  used  in  I,  if  the  neutral  mode  is  subsonic, 
and  satisfies  the  impermeability  boundary  condition  on  r|  =  0,  then  to 
avoid  any  contradiction  the  r i gh t -hand- s i de  must  be  zero  at  the  critical 
point  where  W0  =  c  namely  i)  =  r| j ;  this  requires  (recasting  the  equation 
in  terms  of  the  original  variables) 


Xr 


(4.27) 
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d 

3n 


W0T1 


to  [1+a.c2  +  cm  [* 


+  _ 

a2  (i+xc^+Cn)^' 


=  o  . 

T1  =  T|j 


(4.28) 


This  result  represents  a  further  generalisation  of  the  so-called 
"doubly  generalised  inflexion  condition"  (or  a  "triply  generali-sed  inflexion 
condition"!);  setting  n  =  X  =  0  in  (4.28)  retrieves  the  corresponding 
result  found  in  !.  Note  that  in  the  general  result  for  any  shape  of  cone 
surface  described  by 


r  =  1  +  Xf(C)  .  (4.29) 

merely  requires  the  "X£2"  in  (4.28)  to  be  replaced  by  "X  f(C)" • 

i 

* 

S.  Numerical  RciuIls 

There  are  clearly  many  choices  of  parameters  to  be  made  in  this  study. 
The  strategy  here  will  be  to  carry  out  a  detailed  study  for  one  choice  of 
Prandtl  number  (0.72),  ratio  of  specific  heats  (1.4),  cone  angle  (X  =  1) 
and  Mach  number  (Mao  =  3.8)  for  insulated  wall  conditions;  results  are 
shown  in  Figs.  4-11.  Results  for  cj  are  presented  for  a  range  of  values 
of  n,  at  fixed  values  of  C-  Specifically  results  for  C  =  0-0!  (Fig-  4), 

C  =  0.05  (Fig.  5),  C  =  0.1  (Fig.  6).  C  =  0.2  (Fig.  7),  C  =  1  (Fig-  8). 

c  =  5  (Fig.  9),  C  =  20  (Fig.  10)  and  C  =  75  (Fig-  11).  In  these  figures 

we  just  concentrate  on  growing/neutral  modes.  All  results  may  be  regarded 
as  being  independent  of  numerical  grid.  In  a  number  of  cases  the 
distributions  are  shown  on  two  figures  -  this  is  to  increase  the  resolution 
of  certain  features  of  the  distribution  in  regions  of  small  growth  rates. 

In  all  these  figures,  the  n  =  0  results  are  delinerated  by  a  solidus, 
the  n  =  1  results  by  -  -  -  n  =  2  results  by  —  -  —  -  — , 

n  =  3  results  by  —  --  —  --  —  and  n  =  4  results  by  . 

Neutral  points  are  also  marked  on  the  axes. 
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The  corresponding  planar  results  arc  shown  in  Mack  (1965,  1987a)  and 
I;  the  important  feature  in  this  ease  is  the  existence  of  two  primary 
instability  modes  (and  others,  but  of  considerably  smaller  growth  rate), 
with  the  largest  growth  rate  being  associated  with  the  so  called  "second 
mode" . 

Figure  4  shows  distributions  at  the  location  £  =  0.01,  for 
n  =  0,1, 2, 3, 4.  The  axi symmet ric  (n  =  0)  results  are  very  similar  to 
Mack's  planar  results  described  above  (and  the  axisymmetric  results  for 
the  circular  cylinder  case  close  to  the  leading  edge,  as  considered  in  I). 
The  first  mode  extends  from  close  to  a  =  0  (where  cr  is  approximately 
1-1/M.o  up  to  a  =  0.14  (terminating  at  a  doubly  generalised  inflexion 
point);  we  shall  refer  to  modes  of  this  general  type  as  mode,  I.  A  second 
mode  originates  at  a  =  0.24  (where  cr  =  1,  and  hence  may  be  regarded  as 
a  generalised  inflexional  neutral  mode,  with  the  critical  layer  occuring  in 
the  freestream)  and  terminates  at  around  a  s  0.4  (which  corresponds  to 
a  second  doubly  generalised  inflexional  neutral  mode);  we  refer  to  modes 
of  this  general  type  as  mode  II. 

In  some  ways  the  corresponding  n  =  1  results  for  Cj  (also  shown  in 
Fig.  4)  are  very  similar  to  the  axisymmetric  case.  However  a  third  mode  is 
seen  to  develop,  not  present  in  the  corresponding  axisymmetric  (and  indeed 
planar)  results.  One  important  di st inct ion  between  the  n  =  0  results  and 
those  for  n  *  0  emerges  in  the  limit  as  a  — »  0,  for  which  Cj  f-*  0  if 
n  ^  0  (the  limit  as  a  — »  0,  £  =  0(1)  is  considered  in  Appendix  A), 

although  of  course  the  temporal  growth/decay  rate  ocj  is  nonetheless  zero 
at  a  =  0.  We  shall  refer  to  this  additional  mode  as  mode  1^.  As  a 
increases,  mode  1^  rapidly  disappears,  terminating  at  a  supersonic  neutral 
point  (i.e.  where  cr  <  1-1/VU).  As  a  increases,  a  further  neutral  mode 
soon  emerges,  at  a  second  (supersonic)  neutral  point;  this  mode  is  of  the 
class  I  type  described  previously.  Thereafter,  as  a  is  increased 
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further,  the  n  =  1  distribution  closely  resembles  the  n  =  0  results. 

The  n  =  2  distribution  is  qualitatively  similar  to  those  of  n  =  1; 
however  when  n  =  3,  modes  IA  and  1  arc  seen  to  amalgamate,  although 
mode  II  remains  quantitatively  similar  to  the  mode  II  results  of  the 
previous  n  values.  In  the  case  of  n  =  4,  mode  I  admits  further 
enhanced  growth  rates  although  it  is  still  mode  II  that  possesses  the 
largest  growth  rate  (a  Cj). 

It  is  not  surprising  that  at  £  =  0.01  results  very  similar  to 
corresponding  planar  results  are  obtained,  since  curvature,  generally,  will 
play  a  minor  role  in  the  physics  here;  indeed  a  crude  examination  of  (4.15) 
suggests  that  as  ^  *  0,  the  corresponding  planar  Rayleigh  equation  is 

attained.  Howevcr.as  a  — *  0  and  £  — »  0  a  non-uniformity  is  present; 
this  aspect  is  taken  up  in  some  detail  in  the  following  section,  where 
further  light  is  shed  on  the  additional  mode  1^. 

Figure  5  shows  results  for  Cj  at  (,  =  0.05.  The  axisymmetric  mode 
exhibits  the  same  general  features  as  the  corresponding  mode  shown  on 
Fig.  4,  although  both  modes  I  and  II  have  diminished  growth  rates.  In 
the  case  of  n  =  1,  it  is  to  be  noted  that  mode  1^  has  (already) 
amalgamated  with  mode  I,  whilst  mode  II  has  significantly  reduced 
growth  rates  compared  to  the  £  =  0.01  value,  although  this  remains  the 

more  dangerous  mode.  In  the  case  of  n  =  2,  the  combined  modes  I  and  1^ 

exhibitan  enhanced  growth  rate  when  compared  with  the  £  =  0.01  station, 
whilst  comparing  corresponding  mode  lis,  we  note  that  the  growth  rate  is 
diminished;  indeed,  for  n  =  2,  the  maximum  growth  rate  of  the  two  modes  is 

comparable.  In  the  case  of  n  =  3.  the  first  mode  again  is  more  dangerous 

compareu  to  the  corresponding  mode  at  £  =  0.01.  and  this  trend  is  repeated 
by  the  n  =  4  results,  although  this  value  of  n  indicates  that  an 
increase  in  n  is  causing  a  decrease  in  growth  rate. 
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Figure  6  shows  results  at  the  ^  =  0. 10  location.  These  confirm  the 
trends  observed  previously,  of  a  less  unstable  mode  II.  whilst  the 
(combined)  mode  I  has  a  maximum  growth  rate  comparable  to  mode  II; 
although  it  appears  that  the  n  =  0  mode  II  remains  the  most  unstable. 

Figure  7a  shows  Cj  versus  a  distribution  at  £  =  0.20.  In  this  case 

mode  II  is  barely  visible,  and  so  Fig.  7b  shows  these  modes  on  a  larger 
scale.  Interestingly,  it  seems  to  be  the  n  =  2  mode  I  that  exhibits  the 
largest  growth  rate  and  is  consequently  the  most  important. 

Figure  8a  shows  Cj  distributions  at  £  =  1.  Again,  because  of  the  very 
small  growth  rates,  mode  II  is  shown  on  an  enhanced  scale  in  Fig.  8b. 

Note  that  at  this  axial  location  the  n  =  0  and  n  =  4  (and  above)  cases 
do  not  possess  an  unstable  mode  1  at  all  (whilst  the  n  =  4  mode  II  has 

such  small  growth  rates  that  it  is  not  visible  even  on  the  scale  of 

Fig.  8b).  Overall,  it  is  the  n  =  1  mode  1  which  is  the  most  dangerous. 

Moving  further  downstream  to  -  5  (Figs.  9a,  9b)  we  see  a  "recovery" 

in  the  maximum  growth  rate  of  mode  11.  Indeed,  these  results  show  some 

resemblance  to  the  £  =  0.2  shown  in  Figs.  7a,  7b.  This  is  not  too 
surprising,  given  that,  on  account  of  the  Mangier  transformation 
(Mangier  1948,  Slewarison  1964),  results  as  £  — *  «•  mirror  those 
as  £  — ♦  0  (except  for  a  multiplicative  factor  of  \fT~  in  a). 

This  trend  is  confined  in  Fig.  10  for  £  =  20,  which  may  be 
compared  directly  to  the  £  =  0.05  results  of  Fig.  5.  Notice,  in 
particular,  the  re-emergence  of  mode  for  n  =  1. 

Figure  11  shows  results  at  the  furthest  downstream  location  studied, 
namely  (,  =  75.  The  n  =  0  modes  now  correspond  closely  with  the  planar 
results  of  Mack  (1984,  1987a  for  example),  (with  the  factor  \Z~T~ 
multiplying  a),  whilst  mode  1^  is  visible  for  n  =  1  and  2,  and  the 
union  of  modes  1^  and  I  is  clearly  seen  in  the  case  of  n  =  3. 
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Thus  10  summarise,  these  results,  which  will  guide  us  in  certain 
asymptotic  aspects  in  the  following  section,  we  observe  the  following 
general  features  of  the  stability  of  the  flow:  (i)  results  as  £  — *  ~ 
mirror  those  as  £  — »  0,  except  for  a  multiplication  factor  of  '/IT 
(ii)  Lateral  curvature  has  a  strongly  stabilising  influence  on  mode  II, 
together  with  the  axisymmetric  mode  1  (this  is  in  accord  with  the  results 
found  in  I);  (iii)  there  emerges  a  third  mode,  1^,  as  o  — »  0,  with 
£  — *  0  or  £—♦«•,  in  the  case  of  n  *  0. 

In  the  following  section  we  go  on  to  consider  various  asymptotic  limits 
of  the  system  (4.15)  (or  (4.16)),  guided  partly  by  the  observations  made 
above . 

6.  Asymptotic  Resul ts 

In  this  section  we  consider  a  number  of  asymptotic  limits  of  the 
stability  problem,  to  give  us  a  better  understanding  of  the  details  of  the 
numerical  results  described  in  the  previous  section. 

Perhaps  the  most  intriging  feature  of  these  numerical  results  is  the 

emergence  of  an  additional  mode  as  £  — *  0  (or  £  — ►  «•)  with  a  — »  0.  We 

investigate  this  feature  first. 

6.1  C  — +  Q.  «  =  Oft)  or  C  a  =  Q  (C'h 

Since  the  problem  as  posed  is  basically  equivalent  as  £  — »  0  and 
£  — *  «,  we  consider  just  the  former  limit,  and  later  we  describe,  briefly, 
how  the  results  for  the  latter  limit  can  be  simply  inferred. 

As  noted  in  Section  5.  as  C,  — *  0,  (4.15)  is  seen  to  generally  reduce 

to  the  planar  system  treated  by  Mack  (1984,  1987a  for  example).  However  this 
will  no  longer  be  the  case  if  a  =  0(Q,  since  then  the  denominator  on  the 

left -hand-side  of  (4.15)  no  longer  reduces  to  the  planar  result. 

Specifically,  let  us  write  (consistent  with  (4.12)) 


a  =  ;<s  , 


(6.1) 
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where  it  is  assumed  fit  =  0(1)  as  £  — »  0.  The  results  for  mode 

shown  in  the  previous  section,  together  with  other  results  obtained  by  the 

authors  indicate  that  as  £  — »  0,  c  — »  0  also. 

Partly  guided  by  this,  if  q  =  0(1)  we  set 

c  =  ^cj  +  ^2C2  +  C3c3  +  •  (6.2) 

<P  =  <P0(n)  +  C«Pl  Ol)  +  C2<P2(t1)  +  C3(P3(Tl)  +  -  ••.  (6.3) 

W0  =  W00(I1)  +  C  w0iOO  +  ;2wo2(ti)  +  C3wo3(C)  +  ....  (6.4) 

To  =  T00(q)  +  C  T0i(q)  +  C2T'o2(t1)  +  C3T03(ti)  +  ...  , 

(6.5) 

where  W00(t|)  and  T00(q)  represent  the  planar  values  of  the  velocity 
and  temperature  profiles  repectively,  and  Wqi(t|)  and  Tq i ("H )  etc. 
correspond  to  the  poriurbat ions  to  the  basic  flow  caused  through  curvature. 
To  leading  order,  equation  (4.15)  reduces  to 


d 

3q 


woo  ^oq  •  *ooq  90 
T00  [1  +  n2!  •  M»2  *oo 


=  0 


(6.6) 

t 


i  .  e .  W00  <p0t|  •  *ooq  ‘Po  -  ^o 


(6.7) 


where  k0  is  a  constant.  However  we  require  that  ipo(Tl=0)=  0,  whilst 
90  is  restricted  not  to  grow  exponentially  as  q  — *  ».  Consequently  we 
must  have  that 

k0  =  0  ,  (6.8) 

and  hence 

*?o  =  Ao  *oo  ^  ^  •  (6.9) 

where  A0  is  an  arbitrary  constant  (i.e.  the  unknown  amplitude  of  the 
eigensolut ion) . 

Equation  (6.6)  is  not  a  uniformly  valid  approximation  to  (4.15),  for  all 
q;  specifically  a  breakdown  occurs  when  q  =  0(£’*).  We  define 


q  =  1  +  ^q  =  0(1)  . 


(6.10) 


(i.e.  r|  represents  a  scale  comparable  to  the  radius  of  tbe  cone),  and 
on  this  scale  we  expand  9  as  follows 

9  =  <t>0  (H)  +  C  <&i  (T|)  +  ■  •  •  •  (6.11) 

To  leading  order  equation  (4.15)  reduces  to 


d 


+ 


1 

* 


<P 


0 


(6.12) 


The  solution  to  this,  which  matches  on  to  (6.9)  as  f) 


<t>0  =  Kil 


a  (M.2  .  1>*  t\] 


i  a  (M„2  -  l)*) 


i  s 


(6.13) 


where  Kn ( z ] )  is  the  Bessel  function  of  order  n,  argument  zj  (the 
Kn  (zj)  solution  is  chosen  in  preference  to  the  Im  (zj)  solution  in 
order  that  disturbances  are  propagated  along  characteristics  in  the 
downstream  direction  -  see  Ward  1955,  Kluwick  et  al  1984,  Duck  and  Hall 
1989a, b). 

Returning  now  to  the  q  =  0(1)  layer,  the  0(£)  correction  to  9  is 
given  by 


d_ 

dq 


W01  9oq  ‘  C1  9oq  +  woo  9lq  +  W 00  90  .  W00Tj  9] 
Tr)o  1  1  + 

a* 


•  wolq  9o 


‘00  l  *  +  n2J  -  .^2  ^oo~ 


=  0. 


(6.14) 


and  so  integrating  this  equation  once  we  obtain 


Wol  ^oq  '  C1  9oq  +  woo  9lq  +  woo  9o  *  ^ooq  91  *  ^olq  9o 


=  k  1 


Too  (  1  +  n2  ]  -  .SI*2  W002 


(6.15) 


where  kj  is  an  arbitrary  constant.  9]  must  not  be  exponentially 
large  as  q  — »  »,  and  consequently  we  must  have 
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A0  =  k i  1  -  +  n2 


(6.16) 


together  wi th 


-  cj  A0  W00f|(ri  -  0)  -  kj  Too(0)  1  +  ri¬ 
ft?. 


(6.17) 


However  in  order  to  match  correctly  with  (6.13) 


In  -»  »  =  v°n  'n  = 


=  An  i&  (M^2  -  Q*  Kn  *  *  [itt  (M,2  -  Q* 
Kn'  jift  (NU2  -  1)H 

Eliminating  kj,  we  obtain  the  following  result  for  C) 


(6.18) 


c,  =  1  +  i  ft  (MqJ  -  1)*  Kn"  [  itt  (M^2  -  1)* 


Kn\[ift  (Moo2  -  D* 


?oo  (H-O) 


1  +  n2 


woon  (^O)  l^2  -  1  ■  n2] 

a? 


(6.19) 


The  asymptotic  forms  for  this  expression  in  the  limit  of  large  and  small 


a  may  be  found  readily.  Firstly  as  a  — »  »  we  have 

Cj  _♦  - —  Tqq.^P? - _  +  0(a_1) 

W0on  01=0)  (M„2  -  D1 


(6.20) 


Secondly  as  a  — ♦  0 


Cl  -  ^00°  (n^f  +  m2)  for  ■ 

cl  -  naO°  IS)  *0(62  1085)  for  6=1. 


(6.21a) 


(6.21b) 


In  fact  it  is  quite  easy  to  show  that  as  a  — *  0, 


In  ( c j ) 


a2n  (M,2-  Dn  T00  (0)  22 ~ 2n 
Woon  (0)  ■  [ (n- 1  ) !  I2 


(6.21c) 


Equations  (6.21)  are  precisely  the  (real)  values  found  by  Duck  and  Hall 


(1989b)  for  the  downstream  limit  of  a  non • ax i symme t r i c  viscous  mode 
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(taking  into  account  the  different  scalings  used  in  Duck  and  Hall's  paper). 
Consequently  we  expect  that  as  (,  — *  0/°°,  on  a  scale  smal  ler/largcr  than 
that  of  the  cone  radius,  we  expect  this  mode  to  become  predominantly  viscous 
in  nature,  and  to  be  described  by  triple-deck  theory. 

Distributions  of  I m ( c j }  with  a  are  shown  in  Fig.  12  for  n  =  1,2,3 
(with  (6.20)  also  shown).  Unfortunately  (perhaps)  it  is  seen  that 
I m ( c i )  <  0  for  all  a  (confirmed  by  (6.20)).  From  (6.21)  we  also  have 
that  I  m  ( c  i }  =  0(0  as  C,  — >  0.  Unfortunately,  also,  the  0((2)  and  higher 
correction  to  this  mode  would  require  a  large  amount  of  algebra.  However  we 
are  able  to  make  progress,  in  particular  obtain  an  estimate  for 
Im(c  (a  =  0)}  by  considering  instead  the  limit  as  £  — *  0  of  equation 

(A. 2),  pertainent  to  the  a  — »  0  case.  (With  a  =  o(£)  *we  consider  this 

aspect  in  the  following  subsection). 

Finally,  for  this  subsection,  note  that  the  £  — *  °°  results  may  be 
simply  inferred  from  these  (,  — *  0  results,  simply  be  replacing  the  small 
parameter  '0  in  the  various  expansions  by  the  small  parameter  '  IA£'. 
More  subtle  differences  between  the  (,  <  <  1  and  £  >  >  1  solutions  only 
appear  at  higher  orders. 

fu 2 _ ft  -  Q. _ C  .=i..Q _ Lai _ £-=±ji 1 

The  system  (A. 2)  turns  out  to  be  rather  easier  to  analyse  as  £  — i  0 
than  does  the  corresponding  finite  a  (=  a^')  problem.  We  again  utilise 
expansions  (6.2)  -  (6.5)  (although  see  (6.48)  for  the  <p  expansion). 

To  leading  order,  we  have  for  q  =  0(1)  that 

Vo  =  Ao  ^oo^)  •  (6.22) 

where  A0  is  some  (arbitrary)  amplitude  parameter. 

At  the  next  order  we  have  the  following  system 

^oo  *Plq  ‘  Ao  *ooq  C1  +  Ao  *oo^ 


*ooq  *Pl  +  Ao  *01  *ooq  '  Ao  *o!q  *oo  -  ^1  "^oo  > 


(6.23) 
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where  k]  is  a  constant,  and  we  have  utilised  (6.22).  Setting  T]  =  0  in 
(6.23),  assuming  ^l|q-o  =  ®  •  l^cn 

-Aq  Cj  woon(n-0)  =  kj  T00  (T1=0)  .  (6.24) 


The  boundary  conditions  as  q  — *  «=  must  be  compatible  with  (A. 3).  together 
wi th  (6.22) .  De  fining 

r  =  1  +  rj/C  =  0(1)  ,  (6.25) 

then  we  must  have  an  outer  solution  of  the  form 


out  ?  ,-n-l 
o  =  A  f 


(6.26a) 


where  A  =  A0  +  £A]  +  .  (6.26b) 

and  so  we  must  also  have 

I  out  | 


9111 


n 


=  Vor 


f  =  1  =  - (n+1 )AC 


(6.27) 


(6.28) 


Substituting  this  into  (6.23)  yields 

C1  _  n  Tnn  (t)  =  0) 

1  woorI  (n  =  0)  ' 

in  accord  with  ( 6 . 2 i ) .  In  order  to  estimate  complex  values  of  c  we  must 
proceed  to  higher  orders  in 

At  the  next  order  in  £  we  obtain  the  following  equation  governing  q>2 ; 


^oo  V2q  +  2X  W 00  (p0T^  +  2r)  W00  +  W0 2  Voq 

"wo2ri  Vo  +  wol  Vlq  +  W00  Voq  *  C1  Vlq 

+2wol  ^  Voq  +  Wq1  Vo  *  ^olq  Vl  *  2r|  Wq]^  90 


* 2tj  ci  90T|  "  c2  Vor)  +  ^00  Vl  +  H  ^00  Vo 

-C]  90  -  W 00Tj  92  -  2X  W 00j^  90  -  2q  9j  W007^ 

m  ,  T  „2  T  f  Vo  woo  . 

*tr  Woorj  Vo  =  ^2  ^00  +  n  ^00  j  — 7 - 

0 

+  k!  Tol  .  (6.29) 

We  shall  defer  any  consideration  of  this  equation,  and  move  to  the  next 
order  of  C, ,  which  yields 
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Woo  93q  +  2  XW00  <p  i  ^  +  2  xiW00  <p2rj  +  2  Xti  W00  <p0T| 
+  W00  9lr|  *  C1  92x1  '  2Xc[  90^  -  2q  cj  91^ 

-  H2  C)  90I1  -  c2  9ii]  ‘  2 1]  c2  90q  '  c3  9oi) 

+  W00  92  +  X  W00  9o  +  n  woo  9l  -  C1  91 

■  c  ]  h  9o  '  c2  9o  ‘  ^oot)  93  '  2  X  Woor|  9j 
"  M“  ^olT)  9o  '  2tj  92  W00I^  '  Woorj  9l 
+  wol  92ii  +  2X  Wol  9on  +  2xi  Wol‘  if]r) 

+  h2  wol  •Poxi  +  wo2  *PlT)  +  2X1  Wo2  90T1 
+  wo3  90X1  +  wol  91  +  WqI  n  9o  +  Wo2  90 
Wo3xi  9o  '  2r)  Wo2^  90  -  2X  W0j^  90 

-  2Xxi  W00T^  90  -  Wo2-q  9j  -  2x)  W0i^  91 

-  wolxi  92  =  k2  Tol  +  kl  to2  +  k3  Too 


n2  T, 


00 


[9lT0O  ‘  9o  2T0 1  ]  W0Q  *  cj  T00  90  +  W0i  90  T 


00 


dxi 


‘00 


+  xi'  T, 


ol 


9o  Wqo 


(6.30) 


‘00 


Our  main  goal  here  is  now  to  determine  the  leading  order  imaginary 
component  of  the  complex  wavespeed  c  (we  do,  of  course,  already  know  the 
leading  order  real  term).  Now  since  the  above  equations  just  contain  real 
coefficients,  any  imaginaries  must,  of  necessity,  only  arise  at  a  critical 
point ,  where ,  C  -  Wqo ■  S  ince  c  =  0(£),  this  must  occur  when  x)  =  0(£). 
We  therefore  consider  a  thin  layer  relative  to  the  q  =  0(1)  scale,  namely 

n  =  ti/c  =  o(i)  .  (6.3D 


On  this  scale,  the  expansion  for  9  is  expected  to  develop  as 

9  =  ^90  (fl)  +  ;2  9j  (T1)  +  ;3  92  (n)  +  .  (6.32) 

where  the  9j  are  expected  to  be  normalised  in  such  a  way  as  to  be 
generally  0(1)  quantities. 
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It  is  easy  to  show  that 


V*l>  “  A0  ^oon  *1 


(6.33) 


(where  A0  was  introduced  in  (6.22)),  and  also 


*1  (*t)  =  Ai  W00T1  (n=0)  i)  , 


(6.34) 


where  Aj  is  a  constant,  linearly  related  to  A0 ,  and  we  have  used  the 
property  that  W 00r|q  (r|=0)  =  0). 

After  some  algebra,  the  equation  for  ^(q)  may  be  written 
t^ooq  (0)  *1  •  cil  *  wooq  (0)  ^2 

=  k2  +  3n  Cl  A0  W00T1  (0)  +  A0  w00TiTin  (0)  wOOil(0) 

-  [W00n  (0)]2  A0  q2 


woon  (0)  Toonn  (^)  ci  Ao  ’I 


(6.35) 


where  here  ki  is  an  (arbitrary)  constant. 

If  we  take  (as  we  are  quite  at  liberty  io  do)  A0  and  k2  to  be  real 
constants  (this  is  not  essential  for  our  arguments,  but  simplifies  the 
following  argument),  then  we  now  consider  just  $2'  (where  here  and 
elsewhere  a  superscript  i  denotes  an  imaginary  component).  This  quantity 
is  triggered  by  the  well  known  +  i  71  jump  in  the  logarithm  (Mack,  1984  for 
example)  across  the  critical  layer.  Specifically,  here,  this  is  caused  by 
the  q  dependancy  on  the  r i gh t -hand ■ s i de  of  equation  (6.35)  (k2  plays  no 

role  in  this).  If  (6.35)  is  written  symbolically  as 


[wooq  (0)  i)  -  c |  ]  4>2t]  *  ^oot)  (0)  ^2  -  ^  • 


(6.36) 


*2  =  [wooq  (0)  ^  *  ci] 


[W00 q(0)  q  -  cjJ* 


(6.37) 


Evaluating  this  integral,  taking  just  the  imaginaries  together  with  the 


1  imi  t  as  q 


yields 
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4*2  ~  Aq  (0)  n*C  l]  C|  ^ +  C] 


wo°qgq(0) 

twoon(0)]: 


Toonn(0)  Cl 

Too(0)  O^ooq ( ®)] 


=  B1  [Woon(0)  1\  -  c,]  . 


(6.38) 


i  1  Woo  ^(0) 

where  B  =  tiA0  cj  -  +  C|  - — - , 

wOon(0)  [W<»]: 


(0)  c, 


TOo(0)[W0on(0)]' 


(6.39) 


Equation  (6.38)  then  provides  the  lower  boundary  conditions  for  the  system 
(6.29)  and  (6.30). 

Since  (6.29)  contains  just  real  coefficients  (and  talcing  C2  to  be 
real,  an  assumption  that  may  be  justified  a  posteriori),  then  we  must  have 


<?2  =  B‘  W00  01) 


(6.40) 


where  B  is  defined  in  (6.39). 

If  we  take  the  imaginary  part  of  (6.30)  and  allow  t\  — »  ,  then 


^q  lq 


However  we  require  (on  account  of  (A. 3))  that 


i lii  — *  ~  -  •  ("+l)  Vl„ 


and  so 


~  -  (n+1 )  B 
k3‘  =  -  nB1 . 


(6.41) 


(6.41) 


(6.42) 


Setting  q  =  0  in  (6.30)  yields 


'C1  92q'  |q  =  0  ‘  C3 *  A°  Wooq(°)  *  wooq(°)  ^3  |q  =  0 


(6.43) 


=  -  nB  Too(0) 
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However  from  (6.38)  we  have 

VI,  =  o  -  B'  *'oon<«)  • 


(6.44) 


^3'  It,  =  0  =  c»  B‘ 


(6.45) 


Consequently,  (after  substituting  for  c,). 


-  T  (0)“ji 
1  2  00^  ' 

C  T  =  n  y 

J  Hu  /  n  \  1 


[wooq^)] 


1 


n  T  (0)  W  (0) 
oov  '  oot,t,t,v 


Toot,ti(0)  " 


W0) 


[WOOT,] 


[W°>] 


(6.46) 

Using  the  governing  equations  for  the  basic  flow  allows  this  expression  to 
be  simplified  slightly,  to 


n2  T  (0)2  it 
00 

[^001,(0)] 


1  +  2o  MM2(y- 1 )n 


(6.47) 


In  fact  the  expansion  for  cp  in  (6.5)  is  not  quite  complete  as  it 
stands,  since  the  analysis  of  the  f|  =  0(1)  layer  above  indicates  the 
presence  of  logarithmic  terms;  specifically  we  require 

9  =  <p0(TU  +  C«pi (n)  +  C2<P2(ti)  +  ^3(P3(n)  +••• 

+  log  C  [  C2  <P21  (n)  +  C3  931  (n)  +•  •  •  ]  •  (6.48) 

where  92l(Tl)  =  a21  woo(Tl)  .  (6.49) 

wi th  A21  a  constant . 

A  comparison  of  the  fully  numerical  computation  of  Real  j  c  (a  =  0)  J  , 
with  the  asymptotic  formula  (6.28),  as  £  — *  0  is  shown  in  Fig.  13.  The 
agreement  is  seen  to  be  entirely  satisfactory.  Unfortunately  the 
correlation  between  the  computed  Im  j  c  (a  =  0)  J  and  that  obtained 
using  (6.47)  is  much  less  satisfactory.  However,  this  poor  correlation 
is  not  unexpected  for  two  reasons.  Firstly  accurate  computations  of 
Im  |  c  |  in  this  limit  become  exceedingly  difficult,  as  confirmed  by  the 
quite  complex  asymptotic  structure  detailed  above,  with  both  thin 
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(0=0  (£))  and  thick  (r|  =  0  (1/^))  lengthscalcs  emerging. 

Secondly  the  asymptotic  form  for  1m  |  C|  |  is  achieved  very  slowly  as 
C  — >  0,  at  lcasi  in  one  particular  configuration,  where,  with  n  =  1,  the 
imaginary  wavespeed  has  a  leading  order  coefficient  of  approximately 
3.898  x  105  C3-  A  comparison  between  the  numerical  and  asymptotic  results 
is  not  shown  in  this  case. 

In  the  case  of  £  — >  we  may  just  replace  the  small  parameter 

in  the  above,  by  the  small  parameter  ‘ 1/A.£'.  In  the  following  subsection, 
we  consider  the  behaviour  of  mode  1 .  as  C,  — *  0. 

O—C  -»  o*  -ft  =  Ml h 

The  numerical  results  presented  in  the  previous  section  strongly 
suggest  that  for  the  most  part,  as  £  — *  0  mode  1  has  features  very 
similar  to  the  planar  case,  for  all  values  of  n.  However  there  is  one 
important  exception  to  this,  namely  the  behaviour  of  the  lower  neutral  point 
in  this  limit.  The  planar  case  has  that  c  — »  1  -  1/M„  as  a  — *  0, 
corresponding  to  the  so-called  "sonic"  mode.  However  in  the  case  of  our 
numerical  results,  there  is  evidence  of  a  shift  in  this  neutral  point, 
along  the  positive  real  axis,  and  the  neutral  point  becomes  (slightly) 
supersonic,  with  c  <  1  -  1/M„,  as  £  — *  0. 

A  (sensible)  balancing  of  terms  suggests  that  we  might  look  for  a 
solut ion  of  the  form 

c  =  cn  +  C]  + 

9  =  90  +  C<?1  +  •  ■ 

*0  =  *00  +  C  #ol  + ■ • ■ 

T0  =  Tqq  +  C  Tq i  + .  .  (6.50) 

with 

a  =  ^  a,  a  =  0( 1 ) . 


(6.51) 
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Consequently,  using  (6.53), 
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(^oo  ‘  co)  (^iq  +  9o)  ‘  C1  9oq  *  ^oorj  <Pi  *  1  rj  9o 

K°  [  +  T0 ,  +  2M„2  (W00  -  C0)  cj  -  2M.2  (W00  -  c0)  Wol  ] 

n 

r  (w  -  c  )  <p 

=  Kj  x0  +  a2  Tq  oo  °  dq  .  (6.56) 

1  00 

0 

whc re  Kj  i s  a  constant . 

In  order  to  match  correctly  at  the  outer  edge  (q  — *  *»)  with  (4.18), 
we  must  have  that 


T1  K  (q  )  K  I 
'o  n  v  'o'  o 

Kn'  (%) 


whe  re 


oo 


(6.57) 


(6.58) 


and  qQ  =  a  (2cj)*  (6.59) 

Taking  the  0(1)  terms  of  (6.56)  as  q  — »  ~  and  using  (6.57)  yields 
the  following  non-linear  dispersion  relationship  for  c j : 


(6.60) 


The  integral  (6.58)  was  evaluated  numerically,  and  for  the  conditions 
prevailing  in  all  our  numerical  results  it  was  found  I  -  -228. 4-59. 3i  . 
Equation  (6.60)  was  solved  using  Newton  iteration,  and  results  for 
Real  |  ci  |  and  I m  |  c j  |  foT  various  n  are  shown  in  Figs.  13a  and  13b 
respect ively . 

Notice  that  as  a  — *  »,  (6.60)  predicts  that  one  family  of  solutions 


has  the  property  that 


(6.61) 
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which  is  in  agreement  with  the  a  «  1  expansion  for  c  for  the  planar 


case,  namely 


c  -  1  - 1  /Moo  +  +  0(a4)  . 


(6.62) 


Equation  (6.61)  is  shown  as  an  asymptotic  on  Figs.  14a  and  14b.  Note 
that  the  (real)  family  of  cj  which  also  exist  may  be  found  as  an  exact 


solution  to  (6.60),  namely 


C1  = 


2  M_  a* 


(6.63) 


although  the  importance  of  this  mode  is  not  thought  to  be  great.  The 
complex  families  of  cj's  are  seen  to  terminate  at  a  finite  value  of  a, 
corresponding  to  the  (lower)  neutral  point  of  mode  1.  Notice  that  in  all 
cases  because  Real  j  cj  J  <  0  at  the  termination  point,  these  modes 

I 

correspond  to  supersonic  neutral  modes. 

From  the  result  shown  in  Fig.  14b,  we  are  therefore  able  to  offer  an 
estimate  for  the  position  of  the  lower  neutral  point  of  mode  I  as  £  0. 

In  particular,  for  the  freestream  conditions  considered  throughout  this 
paper,  for  n  =  0  this  position  is  given  by  a.~  0.1  ,  for  n  =  1  by 

a  z  0.20  ,  and  for  n  =  2  by  a  ~  0.295  .  Comparing  these  asymptotic 

results  with  figure  4  in  particular  reveals  a  fair  degree  of  agreement. 

In  the  case  of  C,  »  1,  the  above  results  may  be  easily  transposed, 
by  the  replacement  of  ' by  the  corresponding  positions  for 

lower  neutral  point  are  then  a  -  0.1  (X£)’*  for  n  =  0, 
a  -  0.20  (X£)  *  for  n  =  1,  and  a  -  0.295  (X£)"*  for  n  =  2.  These 
results  are  seen  to  agree  quite  well  with  the  £  =  75  results  shown  in 
Fig.  11. 
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Z _ Conclusions 

In  this  paper  we  have  studied  the  supersonic  boundary  layer  flow,  and 
the  inviscid  stability  thereof,  over  a  sharp  cone  with  adiabatic  wall 
conditions.  The  basic  flow  is  seen  to  evolve  from  one  planar  state,  to  a 
second,  as  predicted  by  the  Mangier  (1946)  Transform. 

The  "triply  generalised  inflexion  condition"  is  derived,  this  being  the 
necessary  condition  for  subsonic  neutral  modes,  and  is  a  (second) 
generalisation  of  the  well  known  generalised  inflexion  condition. 

Significantly,  the  numerical  results  point  to  the  occurence  of  a  third 
mode  of  instability,  not  found  in  similar  planar  and  (more  recently) 
axis  ymme  trie  studies.  An  asymptotic  study  of  this  mode  shows  this  mode  to 
be  linked  to  a  viscous  mode  found  by  Duck  and  Hall  (1989b),  a  study  based  on 
t  r i pi e -deck -  theory . 

The  "sonic"  neutral  mode  is  found  to  be  altered  by  curvature  (and  in 
fact  becomes  a  supersonic  neutral  mode  as  revealed  by  the  asymptotic 
analysis  valid  as  ^  ♦  0  and  £—»«>.) 

Significantly,  our  results  show  that  the  so  called  second  mode  (mode  11) 
is  not  always  the  most  unstable,  at  least  in  the  case  of  non- axi symme t r i c 
modes . 
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In  ihc  case  of  axi symmet ric  modes  (and  indeed  of  planar  modes,  as 
considered  by  Lees  and  Lin  19-46,  Mack  1965a,  b,  1984,  1987  for  example), 
as  a  — *  0,  t he  wuvespeed  c  approaches  the  sonic  limit,  i.e. 


1  ±  1  /Moo  • 


(A.  1) 


In  the  case  of  non  -  ax i symme t r i c  modes,  however,  this  is  no  longer  the  case. 
As  clearly  illustrated  in  the  results,  when  n  *  0,  cj  •/-*  0  as  a  — »  0. 
The  explanation  is  as  follows. 

If  we  (simply)  allow  a  — >  0  (assuming  n  *  0),  then  (4.15)  reduces 


d  (w0-c)  [i+^-C2  +  Ch]  [( +  Cn)  9q  +  C9] 


dq  t, 


won0^2  + 


_  n2  i;2  9 


(W0-c) 


(A. 2) 


As  t]  — »  «>,  this  system  clearly  admits  solutions  of  the  form 


9  ~  n 


(A. 3) 


(together  with  9=0  on  t)  =  0),  which  is  completely  compatible  with 
the  outer  solution,  where  T|  =  0(1),  (r|  defined  by  (4.19))  namely  (4.18). 

Equation  (A. 2)  then  represents  a  reduced  problem  as  a  — »  0,  and  indeed 
a  reduced  form  of  the  triply  generalised  inflexion  point  condition  (4.28) 
exists  in  this  limit,  namely 


d  ^on  O+H2  +  Cn)2 


(A. 4) 


The  system  (A. 2)  -  (A. 3)  was  solved  in  a  number  of  cases  (in  an 
identical  manner  to  the  a  =  0(1)  eigenvalue  system)  and  its  correctness 
was  confirmed.  Notice,  however,  since  the  actual  temporal  growth  rate  is 
ac;,  then  this  still  reduces  to  zero  as  a  — ►  0. 


Fig. 2a  Axial  wall  temperature  distributions, 
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Fig. 3a 
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Fig. 4  Di si ribul ions  of  Cj  with  a,  =  3.8,  £  *  O.Ol 


Fig. 7b 


0.12 


Fig. 8a 
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Fig. 9b  Details  as  Fig. 9a,  enlarged  scale 
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